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ABSTRACT 

We investigate deep inelastic lepton scattering from the nucleon within 
a constituent quark picture, in which the internal structure of con- 
stituent quarks is modeled by meson and diquark dressing. In a co- 
variant framework this structure leads to a breakdown of the factoriza- 
tion necessary for a convolution formulation. We perform our analysis 
in time-ordered perturbation theory in the infinite momentum frame 
where factorization of subprocesses is automatic. Numerical results are 
compared with recent data on valence quark distributions in the pro- 
ton. 
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Many attempts have been made to link the parton picture of the nucleon, which 
is evident in high-energy processes such as deep inelastic scattering (DIS), to the 
familiar quark models of hadrons which successfully describe much of the low-energy 
phenomenology. Most calculations of nucleon quark distributions thus far can be 
grouped into two broad categories — those which use various bag or soliton models 
of the nucleon [I]], and those based on the non-relativistic quark model More 
recently, an alternative formulation in terms of relativistic quark-nucleon vertex 
functions has led to some encouraging results. The advantage here is that one 
deals directly with covariant Feynman diagrams; on the other hand, without de- 
composing the vertex functions into explicit nucleon-quark and nucleon-antiquark 
components the connection between the vertex functions and quark and antiquark 
probability distributions in a general reference frame is not obvious. 

In this paper we suggest how one may combine the advantages of some of the 
existing approaches — maintaining a close link with the constituent quark picture 
of the nucleon while at the same time using a formalism which is fully relativistic. 
We adopt the traditional view that constituent quarks (CQs) themselves have a 
non-trivial internal structure, which can in principle be resolved by a high-energy 
probe. In the impulse approximation, DIS from the nucleon can then be viewed as 
a two-step process, in terms of virtual photon-CQ and CQ-nucleon subprocesses. 

The framework which we use is based on time-ordered perturbation theory in 
the infinite momentum frame (IMF) |Q. Since here only on-mass-shell particles 
are involved, factorization of subprocesses is automatic. Therefore one may express 
the momentum distribution of quarks inside a nucleon (g^r) as a convolution of a 
constituent quark momentum distribution in a nucleon (/q/at) and the structure 
function of a constituent quark (gg). Factorization is generally not possible in a 
covariant framework where particles are off-mass-shell. In any frame other than the 
IMF one must include in addition diagrams with particles moving backwards in time 
("Z-graphs"), which will introduce corrections to any convolution formula. 

To describe DIS from a CQ we consider a simple dynamical model, namely the 
Nambu & Jona-Lasinio (NJL) model P, in which gluonic degrees of freedom are 
integrated out and absorbed into an effective point-like quark interaction, charac- 
terized by an effective coupling constant, g. Since the NJL model is constructed to 
approximate QCD in the low-energy and long-wavelength limit, the CQ distributions 
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will be those at a scale relevant to the model itself (Q^ = Qo ~ 1 GeV^). Explicit 
perturbative gluonic degrees of freedom reappear when the calculated distributions 
are evolved to the higher appropriate for the DIS region. This procedure may be 
questionable since one assumes the validity of perturbative QCD in a region where 
higher order corrections are not negligible. On the other hand, the next-to-leading 
order analysis of Ref. [[7| (which uses Q\ ~ (0.4 GeV)^) indicates 10-20% effects at 
intermediate x when two-loop radiative corrections are included. 

In the NJL model the CQ structure function is determined by the forward quark- 
quark {Q-Q) scattering amplitude [Q, which at zeroth-order in g gives rise to the 

process depicted in Fig. 1(a). In this case the CQ structure function is a delta- 

( s) 

function, and the resulting nucleon quark distribution is given by the function fq/j^: 

q}^"'^^\x) = dpi f^Q%{y,PT) S{1 -xq) = dpi fQ%ix,pT) (1) 

where x = Q^/2P ■ q is the Bjorken scaling variable, y = p ■ q/P ■ q is the fraction 
of the nucleon's momentum carried by the constituent quark, and xq = x/y is the 
fraction of the constituent quark's momentum carried by the struck quark. In Eq. (|lD 
the superscript S denotes a quark distribution associated with a spectator diquark 
system in a spin = or S* = 1 state. The label "bare" refers to a constituent 
quark that is not dressed by Q-Q interactions. 

Beyond the trivial order in g, in the pole approximation the t-channel amplitude 
is well approximated (at least in the pseudoscalar sector) by using meson propa- 
gators, and by diquark propagators in the s-channel. This naturally gives rise to 
a dressing of the constituent quark by meson and diquark "clouds". Figs. 1(b) and 
(c), which has the effect of softening the hard, 5-function, "bare" quark distribu- 
tion in Eq.(|l]). The s-channel diquark-exchange process in Fig. 1(c) in fact generates 
an antiquark component of the constituent quark, which in turn provides a simple 
mechanism by which a non-perturbative sea can arise in the nucleon. The dressing 
corrections will be most relevant in the region of small x, and higher order cor- 
rections involving additional ladders of exchanged mesons and diquarks will not be 
relevant for valence quark distributions, on which we focus in the present work. 

Having calculated their structure within a model, we assume that CQs embed- 
ded in a nucleon have some characteristic momentum distribution, which in the 
present approach is determined by the relativistic quark-diquark-nucleon (QDN) 
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vertex function, $q£)7v The Dirac structure and momentum dependence of the 
QDN vertex can in principle be very comphcated, and our treatment at this point 
is phenomenologicaL Following earlier work 0, Q we approximate the vertex func- 
tion by two structures, one for scalar and one for pseudovector diquarks. A more 
elaborate treatment of vertex functions, for example by solving covariant Faddeev 
or Bethe-Salpeter equations [O], is beyond the scope of the present work, but will 



be dealt with elsewhere 



The CQ structure function and the QDN vertex function then determine the 
valence quark distribution in the nucleon: 

q-'^'\x) = ZgT'^'\x) + (x) - ^^^^gf (x) (2) 

where the three terms correspond to the "bare", meson-exchange, and diquark- 
exchange processes in Figs. 1(a), (b) and (c), respectively. At order one could 
also include the process where a diquark is struck, which would require modeling in 
addition the structure function of an off-shell diquark [|ll|]. In the present analysis 



we do not consider these processes explicitly, but rather assume that they can be 
absorbed into the effective QDN vertex in Fig. 1(a) (see below). The consequences 
for the valence quark distributions will be negligible provided the normalization 
constant Z is chosen to ensure overall baryon number conservation. Direct scattering 
from a pion, however, would enter only at 0{g'^). 

For a scalar {S = 0) diquark we consider a vertex oc (^ — ^')/M, which gives: 

+ 2P -p {rriQ - 171%) - AMrnqp ■ p' - ArrinP ■ p' + AP ■ p' p ■ p') (3) 



where mQ(^D-) ^ 0.4(0.6) GeV is the constituent quark (scalar diquark) mass, and 
sqd = {p + p'Y = ijriQ + Vt) + {tt^d + Pt) /(I ~ ^) is the squared mass of the 
Q — D system. To generate sufficiently hard quark distributions with a structure 
oc / P, |lll one would need more singular behavior in p^ of the vertex function 



^q\)n than the simple dipole form used below. 

Any realistic model of the nucleon must also incorporate pseudovector diquarks. 
We model the quark distribution for an S* = 1 diquark with a vertex 757^, which 
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gives: 
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4P ■ p' p ■ p'^ 



The different 5* = and S = 1 vertices, as well as a larger pseudovector diquark 
mass {rriD ^ 0.8 GeV), explicitly break spin-flavor SU(4) symmetry. This not only 
provides a simple explanation of the A — mass difference, but also the large-x 
F2„,/F2p ratio 

Contributions to nucleon quark distributions from DIS off dressed CQs quarks 
can also be formulated as convolutions, similar to Eq.(|l]). This is only possible, 
however, because of the on-mass-shell condition for interacting particles in the time- 
ordered theory. In a covariant formalism the non-trivial Dirac structure of the 
7*-(5 scattering amplitude and the off-shell dependence of the QDN vertex render 
factorization invalid. This can be demonstrated by writing the CQ structure function 
as: 

(1q{xq) = Tr [(^ + mQ) gQ(a;Q,p2)] (5) 

where the operator qq = I qo + ]^ qi + ^ q2 describes DIS from an off-shell CQ. In 
terms of the coefficients go- -2 the CQ structure function can then be written: 

QQi^Q) = 4mQ go + 4m| gi + Ap ■ q q2. (6) 

For a point-like quark, in the scaling limit one has q2 7^ 0, while go,i = 0. On the 
other hand, for a CQ with internal structure all three coefficients can be non-zero. 
Indeed, dressing by either a meson or diquark cloud gives go = C'(l), gi = C(l), 
g2 = 0{l/q^). In terms of qq the quark distribution in the nucleon is: 

mix) = Jd'p Ti[{lAoip')+^aAUp'))qQixQ,p')] (7) 

where the functions ^0,1 represent the off-shell Q-N interaction. Because of the 
scaling properties of go.. .2, the trace in Eq.(|^ will not be proportional to the CQ 
structure function gg in Eq.^). Therefore Eq.(0) can only be written in convolution 
form if the functions ^0,1 are proportional, however this will not be the case either 
unless the on-mass-shell limit is taken for these. This result is analogous to that for 
an off-shell nucleon bound inside a nucleus 0. However, unlike in the nuclear case 
where nucleon off-shell corrections are not overwhelming ^ 13 1, neglect of quark 



off-shell effects in the nucleoli can lead to significant errors. On the other hand, 
because of the on-mass-shell kinematics in the time-ordered approach, the QDN 
interaction factorizes, allowing 6^^^^q^j^^ and 6^^^^j^^ in Eq.(0) to be expressed as 
two-dimensional convolutions as in Eq.(|l|). 

For the meson dressing process in Fig. 1(b) we restrict ourselves to pions. Dressing 
of constituent quarks by higher mass mesons (e.g. p,u!, etc.) can be considered, but 
contributions from these will be suppressed due to the larger meson masses. The 
contribution to from a constituent quark dressed by a pion is given by: 

5Wgf(x) = f'^ r dp?, fj^%{y,pT) j d'kr qQ^XQ^^PT) (8) 
where /gy^v given by Eqs.(|^) and (0), and 

9Qn \^QQ^{xQ,kT,PT)\^ 



qQ{xQ,kT,PT) 



167r3 XQ{l-XQ){ml - sq^) 
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'k^ + ml{l-XQ) 
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X ( ' ' + xq p| - 2 Pt ■ kT (9) 



represents the 7*-(dressed) CQ interaction. Here ^fg^r = mg/fT, 4 is the Q- 
71 coupling constant p, as given by the Goldberger-Treiman relation, and sg^r = 
{k + k')"^ = {niQ + k^)/xQ + (m^ + (pt — k7-)^)/(l — xg) — p'^. In terms of 
qQ{xQ, kT,PT) the structure function of a (free) constituent quark at rest (Eq.(|])) 
is: qgixg) = dk^ qqixg, kr, 0). 

The presence of transverse momentum pjn in both Jq/n and gg in Eq.(^ means 
that a simple one-dimensional convolution in y alone cannot be valid unless the 
PT-dependence in gg is neglected. This turns out to be a rather bad approximation, 
however, as the pt = results are found to differ by up to 30-40% from the full. 



Pr-dependent calculation of the pion-dressed contributions [0 . 

In Eq. (H) we assume that the integral of gg over kx is regularized by a sharp cut- 
off A, and consequently parametrize the QQn vertex function by a theta-function, 
^QQn = ^ (^T ~ (1 ~ 3;g)A^), with A ^ 0.9 GeV characteristic of a typical chi- 
ral symmetry breaking scale. The contribution of the diquark-exchange process in 
Fig. 1(c) to the nucleon antiquark distribution is similar to the vr contribution in 
Eq.(H), 6'''^^ q^j^\x) 6^^^q^j^\x), but with — > tjid and ^fg^r — >■ gqu ~ 4 taken 
from the NJL model. 

The only additional input necessary for the numerical evaluation of qj^^^^\x) 

( s) 

is the QDN vertex function $gf)^. For this we take a simple dipole ansatz. 



^^QLi^^PT) oc ((A| - ml) I (A| - t)) , where t = -{v\^ x{ml - (1 - x)^)) / 
(1 — x). The cut-offs Ks are chosen to reffect the presumably smaller radius of 
a scalar diquark compared with an S=l diquark: Aq = 1.5 GeV and Ai = 1.2 
GeV. For comparison, we also consider an sqd- dependent dipole vertex function, 
^QDN {sqd + A|)"2, and an exponential type, ^qoat oc exp ((M^ - sqd)/AI). 
These choices are motivated by the resulting symmetry of quark and diquark prob- 
ability distributions under the interchange y ^ 1 — y and rriQ ^ rriD, which one 
would need to satisfy if direct scattering from virtual diquarks was included. Since, 
as mentioned above, these contributions are absorbed into the QDN vertex function 
in Fig. 1(a), this symmetry need not be explicit, and t-dependent vertex functions 
may be used. The main difference between these forms arises at small x, where the 
-SgD-dependent functions give somewhat smaller distributions. The reason for this 
is the 1/x factor in sqd, which at small x serves to suppress the quark distributions, 
which themselves depend on inverse powers of sqd- A full investigation of the effects 
of various vertex functions will be made elsewhere . 

In Figs. 2 and 3 we plot the resulting valence quark distribution in the proton, 
x{uv + dy) = 3x(g™'*-°'' + q^^^^^)/2, as well as the dv/uy ratio. The calculated 
distribution (dashed curve) has been evolved from Ql ^ (0.32 GeV)^ to = 5 
GeV^ (solid curve) using a contour integration procedure similar to that outlined 
in Ref.[0 to invert the Mellin transform. The shaded region represents the range of 
several parametrizations of world data ||15[. Clearly the agreement with the data is 
very good. Even better fits to the data could be obtained by modifying some of the 
parameters in the analysis, such as diquark masses, form factor cut-offs, or the input 
scale Qq. Our choices, however, are motivated by physical arguments and previous 
calculations 0, and we feel the results are encouraging enough to warrant further 
investigation. 

In conclusion, we see that a relativistic constituent quark model formulated in 
the IMF can lead to reasonable valence quark distributions in the nucleon. Inde- 
pendent of any model approximations, the time-ordered approach used here has the 
advantage that factorization of subprocesses is automatic. In a covariant formula- 
tion the off-shell structure of constituent quarks renders the convolution approach 
inappropriate. The present work can also be extended to model the sea of the nu- 
cleon. Already the diquark-exchange process in Fig. 1(c) contributes some 20-30% 
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of the nucleon antiquark distribution at ~ 5 GeV^ |14]. Higher order processes 
involving multi-pion and multi-diquark exchanges, with correspondingly heavier in- 
termediate states, will further enhance the x —>■ region. Finally, the spin structure 



It] of relativistic constituent quarks can be studied within the same framework; the 



details of this will appear in a forthcoming publication [jl 
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Figure captions. 

1. Deep inelastic scattering from (a) a constituent quark in the nucleon; (b) a 
constituent quark dressed by pions; (c) a constituent quark dressed by diquarks. 
The photon four-momentum is g, and P, p and p' are momenta of the nucleon, 
constituent quark and spectator diquark, respectively. 

2. Calculated valence x{uv + dy) quark distribution in the proton (dashed), and 
evolved to Q"^ = 5 GeV^ (solid). The shaded region envelopes the parametrizations 
of world data from Ref. ||15||. 

3. Ratio of proton dy to uy quark distributions. The curves are as in Fig.2. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9406037vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9406037vl 



